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Abstract 

In recent years several techniques have been developed to measure and monitor the pollution of the air. Among 
these techniques, remote sensing using optical methods stands out due to several advantages for air quality 
control applications. A Passive Differential Optical Absorption Spectroscopy system that uses the ultraviolet region 
from 200 to 355 nm of the solar radiation is presented. The developed system is portable; therefore it is practical 
for real time and in situ measurements. The enhanced wavelength range of the system is intended to detect the 
ultraviolet light penetration in the Mexican Valley considering the solar zenith angle and the altitude. The system 
was applied to retrieve atmospheric S0 2 columns emitted either by anthropogenic (power plant) or natural sources 
(volcano), reaching a detection limit of about 1 ppm. The measurement of the penetrating solar radiation on the 
earth surface at the UVC range is presented and the possibility to measure pollution traces of some contaminants 
as 0 3 , N0 2 and aromatic compounds in real time and in situ in the ultraviolet region is discussed. 
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Introduction 

The production of thousands of chemicals has contribu- 
ted to industrial and economic development in many 
parts of the world. This trend however has been asso- 
ciated with the release of new chemicals and possibly 
toxic substances into the environment and food chain 
and adversely affecting human health in many instances. 
The pollutants associated with the anthropogenic activi- 
ties could be inorganic as well organic compounds [1]. 
The air pollution is very complex dynamic phenomena. 
It could be associated with the industrial activities as 
well as with the concentration of the population in large 
towns and cities [2], 

The power plants are using resources like fuel and water 
to provide electricity that is one of the essential needs for 
sustainable development and life. This activity produ- 
ces and discharges all different kinds of pollutants such 
as, gaseous, liquid, electromagnetic fields, and noise which 
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endanger our lives and environments. The amounts of po- 
llutants discharged by power plants, especially air pollution 
are more than assimilation capacity of nature because now 
it is clear that sustaining and assimilative capacity of the 
biosphere though tremendous, is after all finite. The first 
step is to understand the magnitude of emissions from 
each source [3], 

The Differential Optical Absorption Spectroscopy 
(DOAS) has been successfully used for air pollution 
detection in several research works [4]. Important 
advantages over other techniques like gas chromatog- 
raphy, chemiluminescence, and electrochemical, gravi- 
metric, and matrix isolation, are the lower cost and the 
time saving, especially when the sunlight is used as the 
source (passive DOAS). This technique is particularly 
interesting in the monitoring of pollution emissions by 
refineries or thermoelectric power stations. 

Due to the absorption by the 0 3 layer above the earth in 
the UVC interval (200 to 290 nm), the DOAS technique 
has been based on the UVB-VIS range absorption mea- 
surements. It has been used to detect chemicals such as 
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CO, S0 2 , 0 3 , N0 2 HN0 2 , HCHO, CH 3 , CHO, benzene, 
toluene, p- and m-Xylene [5]. 

The absorption lines of these pollutants interfere with 
eachother making more laborious the processing of the 
measured spectra. Usually, the measurement of S0 2 
columns by the passive DOAS technique is done using 
the 290 to 330 nm wavelength range, which is well- 
known penetrates at sea level. In this region the S0 2 
absorption lines interfere with the absorption lines of 
less common species as HCHO, BrO, 0 3 and CIO [4]. 

In the wavelength range from 230 to 260 nm the 
ozone is the principal absorber. Benzene and toluene 
have cross sections one order of magnitude lower than 
0 3 [4]. Ozone atmospheric measurements have been 
previously reported using active DOAS in the range 
from 265 to 345 nm [6,7]. Also, atmospheric 0 3 mea- 
surements using passive DOAS have been done in the 
range from 438.5 to 540 nm. However due to the over- 
lapping of its absorption spectrum with the N0 2 absorp- 
tion spectrum the concentration of both elements is 
measured [8]. 

It is very well documented that the solar radiation in 
the 190-280 nm range is strongly absorbed by the 
stratospheric ozone. Hitherto every report states that 
almost all incoming solar UVC and 90% of UVB are 
absorbed by stratospheric ozone [9]. As a result, there 
are no published experiments measuring the absorption 
levels in the UVC range. However, if there is any irradi- 
ance left, once the radiation passes to the troposphere 
the 0 3 concentration decreases and it would be possible 
to detect the remaining radiation in UVC with a suffi- 
ciently sensitive detector. 

The solar irradiance at the earths surface varies greatly 
depending on factors such as latitude, time of the day, 
month of the year, cloud cover, and haze (aerosols). The 
spectral irradiance in the UVC range on the earth sur- 
face also depends on the ozone column, on the altitude 
above sea level and the solar zenith angle. Very sensitive 
measurements from 280 to 300 nm at sea level, 33 
North latitude and 25.14° Solar Zenith Angle (SZA), 
using PMT detectors, reported the detection of solar 
spectral irradiance of 3xl0~ 4 Wm" 2 nm" 1 [10]. Similar 
measurements have reported a value of two orders of 
magnitude higher at 2500 m altitude, 22° North latitude 
and 7.8° SZA [11]. 

The Mexican Valley is a high plateau in central Mexico 
with a minimum altitude of 2200 m and latitudes bet- 
ween 19° and 20° North Latitude. Under such condi- 
tions the solar zenith angle is almost 0° at noon. By 
what, these values propitiate a higher probability to de- 
tect some irradiance at UVC wavelength range with very 
sensitive detectors as photomultiplier tubes used in 
spectroradiometers or charged coupled devices (CCD) 
used in minispectrometers [12]. In Mexico City the 0 3 



column amounts between 240 and 320 DU [13], if there 
is a possibility to detect some spectral irradiance in the 
UVC region it could be possible to detect traces absor- 
bing in this wavelength interval. So, in order to prove 
the usefulness of the passive DOAS technique in the 
range from 200 to 350 nm in the Mexican Valley, we 
have developed a portable system to measure pollution 
slant columns in this spectral range of several pollutants 
as S0 2> 0 3 , N0 2 , aromatic compounds, etc. 

In this issue the S0 2 emissions around the Tula's Indus- 
trial Complex in the Mexican Valley and around the 
Popocatepetl Volcano are reported. These calculations 
were performed using the developed system. Some discus- 
sion is presented regarding the observed well- structured 
absorption spectra in the range of 200 to 290 nm. 

Materials and methods 

The system consists of a fused silica lenses telescope 
84-UV-25 (Ocean Optics) with a 25.4 mm diameter, 
focal length of 100 mm. This telescope is especially sui- 
table for light collimation from a distant source. The 
collimated light is guided to the spectrometer using an 
optical fiber (QP600-2-SR, Ocean Optics) with a 600 m 
diameter core and 2 m long. The spectra measurement 
was performed using a miniature spectrometer (USB 
2000, Ocean Optics) with a diffraction grating of 2400 
lines/mm. A cylindrical lens is incorporated to the de- 
tector to increase light collection efficiency focusing a 
1 mm high slit onto a considerably shorter (200 \im) de- 
tector elements. The CCD detector (ILX511A, Sony) is 
coated with an UV anti- reflecting film to reduce the re- 
flection losses and therefore increasing their sensitivity. 
The spectral range of the spectrometer is in the 190 to 
355 nm range with a resolution of 0.5 nm. This system 
shows a high signal to dark current random noise ratio 
in the 190 to 290 nm range (over 3 times for measure- 
ments under pollution measurements and about 12 
times under clean sky measurements). The system was 
checked and no overlapping of different diffraction 
orders and the scattering of the incident light which 
could affect the readings in the UVC region was 
detected. To reduce the background noise of the mea- 
surement an averaging of 3 to 5 spectra was performed. 

The coordinates of the measuring spots were deter- 
mined with a GPS locator (GT 37231). All the information 
inputs a laptop were the developed software processes the 
data. The coordinates given by the GPS locator were 
superimposed on the Google satellite images using the 
Google-Earth-KH-80.11b libraries to show the actual mea- 
suring points on a satellite photograph. 

The attenuation of the irradiance I 0 (X) is described 
by the Beer-Lamberts Law 1(A) = I 0 (A)e- a(X) 1 c(s)ds . The 
attenuation depends of the concentration c(s) of the 
different compounds of the light path media [14]. The 
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absorption cross section o(X) defines the light absorption 
to a specific wavelength. The dependence of the cross 
sections <r(A) on the temperature and pressure has not 
been considered in this work. The gases columns are 
evaluated by using the mixing ratio method [15]. The 
light absorption cross section of a specific pollutant can 
be obtained from reported data on literature or from la- 
boratory measurements. 

Usually, it is very difficult to obtain separately the pol- 
lutant concentrations c(s) and the optical length /= J ds. 
Therefore, these parameters are combined to form what 
is called slant column density S= J c(s)ds expressed in 
molecules/cm 2 [16,17]. 

Actually, the measured absorption involves all the present 
molecules that absorb in that specific wavelength. The ab- 
sorption of several substances i will be then described by 

the following expression: I(X) — Io(X)e ? ( } , where <X;(A), 
Ci(S) and S t are their absorption cross section, its specific 
concentration and its slant column density, respectively. 
The sunlight spectrum that reaches the Earth surface shows 
slow and fast variations. The slow variations of the 
spectrum are determined by the light source and by the 
scattering. The fast variations are determined by the ab- 
sorption of molecules. The spectrum of the cross sections 
can be divided in a slow variations of (A) and fast variations 
parts of (A). Therefore, the absorption cross section for each 
component will be a^X) = of (A) + of (A). 

The slow variations part of the measured light intensity, 
which depends on Iq{X), can be approximated by a polyno- 
mial. Using some polynomial adjustment method or by 
using high-pass digital filters, for instance the Savitzky- 
Golay [18] a solution for fast variations part of (A) could be 
obtained [4,19]. 

For the detection of only one gas component in a 
spectral interval, if the scattering is small compared to 
the slow variations of radiation I 0 (X) the pollutant slant 
column density can be expressed as: 

o F (X) o F (X) 

Where D(X) is the Optical Density of the medium 
[4,15]. The cross sections of the gas trace at the atmo- 
sphere could be considered independent from their 
placement; this means that they are independent from 
their altitude [20]. 

High resolution absorption spectrum for the S0 2 has 
been reported [21]. The reference spectrum for the 
evaluation of the measurements is needed and it must 
consider only the fast variations of (A) and the spectral 
deformations produced by the spectrometer and deter- 
mined by the spectrometer transfer function (STF). In 
practice, to obtain the STF the spectral lines emitted by 



a Hg lamp could be used because, these are narrower 
(-10 pm) [22] than the typical resolution of the used 
spectrometer (0.1-1 nm). The measured spectrum is a 
good approximation of the spectrometer transference 
function [17]. Mathematically the measured spectrum 
/ (A) can be expressed as the convolution of the Hg lamp 
emission spectrum l(X) with the spectrometer transfer- 
ence function F(X) = I{X') * H sp (X'), where H sp (X) is the 
spectrometer transfer function. The STF was measured at 
334.1 nm because at this wavelength the "cleanest" signal is 
obtained. 

DOASIS software, developed at Heidelberg University 
[23], was used to perform the convolution of the trans- 
ference function of our instrument. The high resolution 
S0 2 spectrum is used to obtain our reference spectrum. 
To observe if there is any shift of the S0 2 reference 
spectrum related to the convolution process, the S0 2 
absorption spectrum with a Xe lamp (Newport, model 
969607) was measured. A constant shift of 0.98 nm was 
obtained and it must be considered when the trace gas 
concentration calculation is performed. Figure 1 shows 
the S02 reference and sample spectra. 

A 2.5 mbar, 3 cm long, S02 cell was used to validate 
our DOAS system. We used a Xe-lamp as radiation 
source. The measured column agrees within a +- 4% 
error. To determine the density of a measured column a 
non-linear adjustment is done by using the Levenberg- 
Marquardt method [24] using the measured I 0 (X) and 
the reference spectra Io(X). A dark spectrum, previously 
captured, is subtracted from I(X) and I 0 (X). A residual Ai// 
is obtained from the fitting process, in our case in the 
order of 2.5xl0 17 particles/cm 2 , which can be used to 
calculate the detection limit of our system trough the 
relation S min = ^ [17]; hence the detection limit is the 
order of 1 ppm. 

The spectrometer is connected to a laptop. A program 
was developed in Labview to capture and filter the mea- 
sured spectra to perform the non-linear adjustment and 
to identify the amount of polluting gas. This software 
consists of several programming stages. First of all, a 
driver for the Windows operative system is instal- 
led from the spectrometer that comes from the 
OOIDDrv32.dll (Ocean Optics, 2011) library. Once the 
hardware is recognized a Sarvitzky-Golay filter program 
filters the reference I 0 (X) and measured I(X) spectra. A 
dark spectrum is subtracted from and I(X) to reduce the 
thermal noise of the instrument itself. Immediately, the 
relation of these spectra is calculated and is applied the 
logarithm function. Spectra are geographically refe- 
renced throughout the journey to integrate all the S0 2 
measured columns. A mouse-type GPS was used with 
the NMEA protocol to link it with the program, to do 
so, an API for Windows is used which is provided by 
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Figure 1 Reference and Sample spectra to determine the S0 2 content. 
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Labview under the VISA name, which allows to read the 
data from the GPS ([25]; Visa, [26]). 

When the DOAS technique is used, the path length of 
light depends on the number of scattering events that 
take place until the photon reaches the spectrometer. 
Some methods to obtain at least the statistics of photon 
propagation are based on the simulation Monte Carlo 
method [17] to be able to validate statistically what ha- 
ppens with a large amount of photons it is necessary to 
simulate such events unfortunately this process requires 
a computing robust system and the simulation process is 
not performed in real time. That is not practical for 
insitu measurements. Taking this into consideration, in 
this work the column densities are indicated in ppm*m, 



using mixing ratios [27,28]. The S0 2 flow is calculated 
in tones by day [14], using column densities measured 
by the system, which must be multiplied by the distance 
of the traverse route using the following formula: 
ton/day = 0.0OO23j%SCh(ppm)*X*VsenddX Here X is 
the traverse distance in meters from the point Xj to x w V 
is the wind speed in m/s, 0 is the angle between vectors 
X and V. For the real time in situ measurements around 
the industrial zone and the volcano, the telescope was 
fixed on a car for the traverse. Five spectra were ave- 
raged in each data entry to minimize the noise effects. 
The spectra are taken every 3 seconds. The system takes 
firstly the dark spectrum to reduce the noise produced 
by the spectrometer and the skylight spectrum I s /dW is 
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Figure 3 Optical density observed from 200 to 350 nm. 



measured. Finally, the sample spectra are taken during 
the traverse around the industrial zone. 

Results 

The reference and measured spectra used for the deter- 
mination of S0 2 content in the 280-320 nm wavelength 
range are shown in Figure 1. The reference spectrum 
was obtained according to the method described earlier. 
The absorption bands of S0 2 are very well defined in 
this interval. 

The optical spectra in the 200-350 nm wavelength inter- 
val with and without contamination are presented in 
Figure 2A. The skylight spectrum is detected with a "clear" 
sky, which means, without contamination plume. It is 
clearly observed the attenuation of the spectral irradiance 
on the Earth surface due to the pollutants absorption. 
More than that, it is perfectly observed that a certain "non 
despicable" irradiance is detected in the UVC region. The 
UVC wavelength interval has been enlarged and presented 
in Figure 2B where an absorption band structure is ob- 
servable. The attenuation in the plume in the range bet- 
ween 200 to 290 nm of the measured spectrum I(X) could 
be related to the 0 3 concentration, which is the main 
atmosphere component absorbing in that region. 

The optical density is presented in Figure 3. The cha- 
racteristic S0 2 absorption bands are very well defined in 
the spectrum range between 290 and 320 nm. Besides, a 



band structure related to scattering and absorption in 
the presence of other molecules in the UVC range can 
be observed. Usually, this interval is neglected in most 
studies. This result shows the importance to develop 
some similar algorithms to evaluate the concentration 
of species absorbing in this range, such as tropospheric 
ozone. 

After the high-pass filter is applied and the spectrum 
restricted from 290 to 315 nm range, the differential ab- 
sorption S0 2 spectrum is obtained and it is presented in 
Figure 4. 

The developed system was used to measure both, 
anthropogenic and natural S0 2 emissions. Two sites 
with significant S0 2 concentration were chosen near 
Mexico City. 

Measurements at Tula hidalgo refinery 

In Tula Hidalgo, Mexico, two nearby industries are located: 
"Miguel Hidalgo" PEMEX oil crude refinery and CFE 
thermo power station "Francisco Perez Rios". The satellite 
image is presented in Figure 5A. Once the traverse was 
completed, the S0 2 columns are calculated using a wave- 
length interval from 290-350 nm. The variations in ppm*m 
can be observed in each column and the results are shown 
in Figure 5B. Both industries are responsible for the emis- 
sion of nearly 350 t/d of S0 2 [29]. In this sense, it is a con- 
venient place to test the equipment. The wind direction is 
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Figure 4 S0 2 differential absorption spectrum in the 290 to 320 nm interval after filtering. 
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Figure 5 A. Satellite image from Miguel Hidalgo refinery in Tula. The S0 2 polluting chimney, the wind direction, and the traverse route 
are presented. B. shows S0 2 measured columns with passive DOAS system around Miguel Hidalgo PEMEX refinery in Tula Hidalgo. 



presented on the satellite image, Figure 5A. Also, the 
traverse route is shown on both, the satellite image 
(Figure 5A) and the corresponding S0 2 column distribu- 
tion (Figure 5B). The starting and ending points of the 
traverse are indicated. The result presented in Figure 5 
(A and B) corroborates that the distribution of the S0 2 
columns with higher concentrations are related directly 
with the wind direction. 



Measurements at Popocatepetl volcano 

Popocatepetl volcano is located at the boarding limits 
between Estado de Mexico, Puebla and Morelos. Its 
coordinates are 19.023° North Latitude and 98.627° West 
Longitude. A traverse from San Buenaventura Nealtican 
to San Nicolas de los Ranchos in Puebla was done; it is 
shown in Figure 6A. It has to be noticed that the tra- 
verse direction is from east to west, as the starting and 
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Figure 6 A. Satellite image from Popocatepetl volcano in which the wind direction and the traverse route are depicted. B. shows S0 2 
measured columns at Popocatepetl volcano with passive DOAS system. 



ending points are indicated. The measured S0 2 columns 
are shown in Figure 6B. The first part of the traverse is 
oriented almost radially to the volcano coordinates. At 
the starting point the system was calibrated. As the 
measuring point approaches the volcano the S0 2 co- 
lumns readings increase, as it is expected due to the dif- 
fusion and the dilution of the pollution cloud. 

At some point of the traverse, indicated with an arrow 
on Figure 6A, the direction changes and it becomes 
nearly perpendicular to the volcano's direction. In this 



second part of the traverse the measured values of S0 2 
columns remain practically constant. 

Discussion 

Usually, the measurement of S0 2 columns by the passive 
DOAS technique is done using the 290 to 330 nm wave- 
length range. In the wavelength range from 230 to 
260 nm the ozone is the principal absorber. Albeit the 
arrangement of our spectrometer shows less resolution 
than a similar one with a range from 290 to 355 nm, 
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from Figures 5 and 6 it can be clearly observed that the 
developed system has still enough resolution to determine 
the S0 2 pollutant distribution and, after Figure 3 (for alti- 
tude and latitude conditions as in the Mexican Valley), 
also could be used to simultaneously measure 0 3 columns 
in the 230 to 260 nm range. 

Still is necessary the measurement of the wind direc- 
tion and speed in the plume at the same time with the 
column measurement to estimate the flow of the S0 2 ; 
this is an open problem. If the wind conditions remain 
steady as the columns measurements along an appro- 
priate traverse are performed, a good approximation of 
the emitted concentration and the pollutant flow could 
be calculated. 

Conclusion 

In this work we describe a passive DOAS system deve- 
loped to determine S0 2 polluting columns. Portability of 
the developed system makes it practical for real time 
measurements and it allows the S0 2 columns measure- 
ments in situ. It could be applied to develop some mo- 
dels of pollutants distribution considering the source, 
the diffusion and the wind conditions. The direct meas- 
urement of the pollutant columns provides useful infor- 
mation to evaluate the pollution level in real time. Also, 
DOAS technique avoids contact with the samples pro- 
viding a better interpretation of the S0 2 columns 
measurements. 

Attenuation in the plume of the reference spectrum I 0 
(A) in the wavelength interval from 200 to 290 nm is 
presented in Figure 2B. It should be emphasized that a 
certain level of irradiance in the UVC region is detected 
on the earth surface. It could be related to the altitude 
and latitude of Mexican Valley, as well as with the con- 
tinuous improvement of optical detectors. The absorp- 
tion in the UVC range could be related to the 0 3 
concentration. It leads to a possibility to perform DOAS 
in this wavelength band applied to ozone detection and 
monitoring at Mexico City region, due to its altitude and 
latitude. However, much more research work should be 
done to corroborate this hypothesis. 
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